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ABSTRACTS 


By meQ.fis of aevial photography and MSS-LAWSAT data ^ a land use/ 
land cover classification was applied to the Tubardo River coas'bal 
plain. The following classes were identified: Coal related areas j 
permanently flooded wetlands^ periodically flood/ed wetlands^ 
agricultural lands ^ hare soils, water bodies, urban areas, forestlands. 

1. INTRODUCTION 

The acquisition of land use and land cover* information of a costal 
region is among the first steps toward a rational development 
planning. This information is important because these land covers 
occupy a transitional zone bordering land and sea, acting as a buffer 
to the environmental interactions between these two major divisions 
of the coastal ecosystem (Zenkovich, 1967) . As ecotone communities 
between sharply changing environments, the recovery of the natural 
vegetation after a disturbance (burning, draining, water logging, etc) 
may not happen, and an erosional process takes place until another 
equilibrium point is reached (Tomlins and Thomson, 1980; Godfrey 
and Godfrey, 1974) . The bordering tide affected communities are 
specially important as nursery niches for many of the aquatic fauna 
and as primary producer to the estuarine detritic food chain (Odum 
and Cruz, 1967; Gallagher et al . 1980; I-Jhite et al. 1978). 


* The distinction between land use and land cover is not necessary 
for a remote sensing problems, so they will be referred as layrd 
cover. 
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In despite of the fact that the composition of these communities 
be relatively simple, the determination of their localization and 
areal extent by field methods are not an easy task due to the 
difficult access to some areas with extremely soft soils. Several 
authors report that remote sensing techniques are successfully used 
in coastal land cover mapping (Carter and Schubert, 1974; Gallagher, 
1974; Bartlett et al. 1976), Since it presents a synoptic view of a 
subject area, the identification, localization and aerial extent 
calculation is fast and cost effective with a minimum ground 
checking. Some attempts have been made to extract quantitative 
information such as green biomass on coastal zones, by means of the 
vegetation reflectance characteristics in the red and infrared 
spectral bands (Gallagher et al, 1972; Bartlett and Klemas, 1979). 

In order to evaluate the capacity of the remote sensing data 

available at INPE to detect coastal land cover changes, a research 

project proposed the analysis and comparison of the Tubarao river 

coastal plain’s land cover in 1978 and 1981. An area of approximately 
2 

600 km was chosen as study area and it encompasses the cities of 
Tubarao, Laguna, Jaguaruna (Fig.l). 

The Tubarao river coastal plain is formed by quaternary fluvial 
deposits which started to accumulate when the once Tubarao bay 
was closed by sand spits formed by the sea, along a string of 
relief islands, about 25 km off the shore of the bay (Reitz, 1961; 
Guerra, 1950) . Nowadays, the enclosed bay is mostly occupied by 
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marsh covered organic soils with an average level of 2 in. A few 
remaining water hodies are present on the coastal plain. Marine 
terraces are present in the South of the Garopaba lagoon and a 
fluvial terrace is noticed near the Tubarao river mouth. Granitic 
outcrops divide the coastal plain in an inselberg like way with 
height ranging from 20 m to more than 400 in. 

The recent history of this coastal plain presents many strong 
environmental disturbances, mainly on water quality and availability. 
As consequence, agricultural areas are substituted by natural vege_t 
ation in some places while the opposite may be happening somewhere 
else. Total vegetation deterioration occurs in places of severe water 
quality impact. 

A brief description of the main environmental disturbance on this 
area is presented below: 

~ From 1964 to 1978, pyrite-rich coal refuse was deposited in the 
Northern part of the Estiva marsh field. A large area of 
totally killed and heavily damaged marsh lays along the acid 
drainage. 

- By 1980 this p^nrite-rich spoil began to be explored by a 
pyrite concentration plant. The area of vegetation damaged 
by the acid drainage encroached while a new carbonous refuse 
site was built damaging another place dovmslope. 

- The Camacho lagoon outlet, closed by the time, was artifficia^ 
ly opened to drain a large amount of trapped waters on the 
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Garopaba area when the Tubarao river had a great flood in 
197A. Natural sand deposition completely closed it by 1980. 

- The Tubarao river, which used to meander across the plain, 
was straightened and lowered to prevent such major floods. 
This work was mostly accomplished by 1980. 

- New agricultural areas on former wetlands are rather 
frequently installed, wherever investments are made on 
constructing ditches and drainage channels. Agricultural 
areas are also abandoned because of soil impoverishment, rise 
of the water table and silting of the drainage network. 

As this paper will present the results from the 1978 analysis, its 
objective is to map the land cover of the study area by means of the 
available remote sensing data and to evaluate the results qualitativ^e 

ly. 


2. METHODOLOGY 
2.1 “ Materials 

In the development of this research it v?as used: 

- A set of 52 color infrared (CIR) aerial photographies at 
the scale of 1: AS, 000, with longitudinal overlap of 60% and lateral 
overlap of 30%, obtained in Augus t/September of 1978; 

“ MSS LANDSAT data CCT (path/row annotation 178/32 from 

April 2A, 1978) with digitized information of the spectral radiance 
measured at every 57 x 79 m picture element in four spectral bands: 
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.5 “.6 ym, ,6-. 7 pm, .7-. 8pm, .8-1.1 pm. 

- DF Vasconcelos pocket stereoscope with 2 x magnification 
power ; 

- INPE’s Multispectral Image Analyzer System - Image-100 
(I-lOO) with its ovm hardware for image processing; 

- The following software for digital jjnage processing developed 
at INPE; 

- MEDIA K - An unsupervised K-means classification algorithm 
(Dutra et al, unpublished) 

- MAXVER - A maximum likelihood decision rule based 
supervised classification algorithm. (Velasco et. al., 1979) 


2.2 - Methods 

A) Air photography analysis. 

A classification system based on Anderson et al , 1976 
and Cowardin, 1982, was built in order to hierarchically organize the 
land cover classes observed on aerial photographies. The classes were 
charact eryzed by their tone, texture, context and, specially for man 
made classes, size and shape. Some of the classes were named based on 
the photointerpreter experience while some others were kept in a 
descriptive form until a ground check was made. The classes obtained 
by this procedure will be mentioned as information classes. 

Land cover classes as well as drainage networks, roads ^ 
and main reliefs were outlined on a sheet of poliester for each aerial 
photography. The poliester sheets were set together with the aid of a 
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topographic map, but without the sufficient number of ground control 
points to be considered as a semi-controlled map. Two final products 
were obtained: a map describing the land cover patterns and the maXn 
roads and another with the drainage features and the borders of the 
relief forms. 

B) MSS-LANDSAT data analysis 

The below described procedures were interactively 
evaluated at the I-lOO TV monitor at the scale of 1:104,100. The 
intermediate results were hard copyed as print-outs at the same 
scale of the photo mosaic. 

As many of the classes obtained from the aerial 
photographies were identified from context clues, MSS-LANDSAT data 
classification approach avoided to use the ground truth map as the 
unique information source for the training phase. Spectral information 
was extracted from LANDSAT data themselves using an unsupervised 
classification algor itlim, the MEDIA-K algorithm, which allowed to 
divide the image into 8 spectrally homogeneous classes. Using the 
aerial photography derived map as ground truth, each spectral class 
was analyzed separately. 

As each of the spectral classes represented two or 
more third level information classes, the MAXVER algor itlim was used 
in an attempt to subdivide them. The training samples for the MAXVER 
algoritlim were chosen in places where the information classes presen_t 
ed a spectrally homogeneous classification. The hard-copy of this 
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procedure, vith 18 cXaBses, vas again checked against the ground 
truth raap. It vas noticed tha?. many hierarchically superior* 
information classes could be represented by the addition of two or 
three of the obtained classes. Most of the final classes were thus 
obtained by the addition of the fonuer procedure’s class samples, 
except for the bare soil class, which was kept as three separate 
spectral classes represented under a single color display or print~ 
-out s^aubol, due to its intrinsic spectral differences as consequence 
of soil origin and water status. 

The final MSS— LANDSAT data classification out— put Vfas checked 
against the ground truth map in order to qualitatively evaluate the 
performance of the developed image classifier. 

3. RESULTS AND DISCUSSIOMS 

3.1 - Classification svstem: 

- - — - - - ■ ■ 

The classification system (Table 1) was hierarchically 
built so that it could be applied to both aerial photography and 
MSS-LAMDSAT data (Anderson et al., 1976). A brief comment will be 
made about those classes which need an citplanation to better define 
the classification system. 

The coal related class was established as a first level 
class because of its importance as land cover change agent. The 
detailed subdivision of this class is duo to the differential 
pollution potentials of its categories as well as to the differential 
MSS.LANDSAT identification capability for these categories, according 
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to their spectral characteristics and the presence of water. 


An approach, closer to the classification system, proposed 
by Cowardin, 1982, was applied to the wetland class. The water regime 
of the area was the main criteria for the level II subdivision. 
Fortunately, the MSS-LANDSAT infrared channels work out very accor^ 
ingly to this approach. The Typha dQmu^guot'ioto is frequently found in 
the wettest parts of the coastal plain, except near the Santo Antonio 
lagoon, where the salinity reaches higher concentrations, and 
the Typha dominguens’is is substituted by Spavt i na sp. stands. 

On the sandy marine deposits, along the shore, depressions between 
deposition rills are ioccupied by fresh or brackish water marshes 
named by beach hydrosere. It corresponds to the ’’Etapa paludosa" of 
the Reitz’s description (Reitz, 1961). The coal stressed vegetation 
class presents acid damaged plants ranging from totally killed to 
heavily chlorotic stands (Murtha, 1976). The floating hidrosere is 
composed mostly by Eichornia sp. The pure sedge class, represented 
mostly by Sayx>pus giganteoiw stands can cover large periodically 
flooded areas. However, it is present on drier soils which are 
covered by pastures, after the original stand is burned. The under 
MSS-LANDSAT resolution mosaic, formed by irregular fire and the 
fast recovery of S, gigajitcouB was categorized as sedge plus pasture 
class. The Estiva march class is originally height dominated by an 
non- identified Compos it ae. Frequent slash and burn turns it into 
a grass covered meadow which is quickly substituted by a Cyperacea> 
Juncaceae community. 
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The dreirted pasture lands, sometimes with cultivated species 
such as Bvaohiaria op and S^itar*i<x op, have their water regime 
artificially controlled. The frequently flooded pasture lands occupy 
once B, gigant 60 U 8 dominated areas and arc periodically burned to 
avoide the regrowth of the original community. The meadows, periodically 
flooded, are composed by a very diverse herb community. Burning is not 
frequent on this class. 

The "jundu” class is the fb-ed dunes vegetation dominated, at 
the study area, by the fire-resistant Butia capitata and a large variety 
of shrubs, mostly Myrtaceae. The moving dunes communities are stands of 
xeromorphic herbs and grasses occupying the between-dune spaces or even 
pioneering the dunes themselves, 

3.2- Ground truth map 


The ground truth map, obtained from the air photographies 
analysis, is presented on Figure 2. Land cover patches, smaller than 
four MSS-LANDSAT picture elements, were not represented unless their 
contrast against the surroundings was very ligh. 

3.3 - MSS-LANDSAT derived map 


The final product of the MSS-LANDSAT data classification presented 
eight classes, i.e., coal related areas, permanently flooded wetlands, 
periodically flooded wetlands, vegetated agricultural lands, bare soils, 
water bodies, urban areas and forestlands. An alphanumerical representation 
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of these classes can be seen at Figure 3. The performance of each 
class is given below: 

. Coal related areas: all but one of the level III classes 
were identified on the MSS-LANDSAT map, correctly positioned and 
with borders very similar to the actual. The ash slurry ponds were 
onimited for the permanently flooded wetland class. Shallow waters 
were commissioned into the coal related class and this is observed 
along the shores of the lagoons. 

. Permanently flooded wetlands: a few commissions of relief 
shadowed areas as well as the ash slurry ponds were only erros 
revealed by this class. 

. Periodically flooded wetlands: all of the pure stands 
of S. giganteous were ommited from this class to the vegetated 
agricultural lands. As one of the main MSS- LANDSAT features used 
to discriminate the wetland classes was the antagonic response of 
green biomass and water on the infrared channels, the complete 
cover provided by the dense, 2-3m tall S. giganteous canopy, hide 
the water and moves this class’ spectral signature into drier 
vegetated classes decision space. The other two classes were 
satisfactory. This class also presented full commission of the 
frequently flooded pas tur elands , 

. Vegetated agricultural areas: unf ornately the discrimination 
of the main agricultural land uses of the study area, rice fields and 
pasturelands was very difficult for the utilized date. In April the 
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rice fields ere harvected and thv green parts are being grazed by 
cattles. After several failed attempts to separate these classes 
with the interactive algorithm IIAXVER, they were considered as a 
single class which included the others class (corn, cassava, vegetables, 
etc). Visually this class has the lowest ommission error, but 
commissioned areas of microtextured mixed unvegetated and vegetated 
classes and all the jundu area. 


. Bare soils: the multimodal distribution of this class’ 
intensity histogram for any of the MSS-IANDSAT channels, makes 
necessary the subdivision of the class. The color of bare soil 
fields on a false color infrared composite simulation on the I-lOO 
TV monitor ranged from black, dark green, ..., green, light 

green and white. The white fields were associated with the marine 
formations near the city of Jaguaruna, The others are known as peat 
soils with varying organic matter content and were divided arbitrarilly 
into two classes. They were added into a single class for color 
representation of the I-lOu system (maximum of eight colors) . A 
classical confusion between this class and urban areas class is 
present and can be seen in the South of Tubarao city. Dark soils 
were also omitted for wetland classes. 


. Water: Due to the characteristically dark tonality of the 
water response on the MSS-LANDSAT infrared channels in contrast to the 
brighter land classes, the mapping of the water bodies was sucessful 
for all lagoons larger than 5 ha. The tide influenced lagoon shores 
were frequently omitted for coal related class. 
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.Urban areas: the larger urban scgiuents vere reasonably 
mapped, but the urban fringe areas as well as the settleincnts along 
the main roads were mostly omitted for the soil and vegetated 
agricultural lands classes, respectively. 

. Forestlands: The absence of large stands of forests in 
the study area hampered the image sampling and good spectral 
characterization of this class. The results identified only those 
forest stands on sunny slopes of the rock outcrops. The partial 
soil cover of the "jundu*' conmiunity is the reason for its 
classification as agricultural lands, 

/u CONCLUSIONS 

The MSS-LANDSAT derived classes presented varying identification 
performance and, of course, varying information reliability. For 
non-vegetated classes, coal related areas and waterbodies are, 
except for specific errors, accurately described. The sand dunes 
and bench area can be described by the unclassified stripe, bordering 
tlie ocean in the print-out. One can rely only on the main cities 
cores description by the urban areas class. The same applies to tbo 
bare soil class which worked out well only on large ploughted 
fields. The worst performance was noted on the forestlands class. 

The other three vegetated classes have a fairly good classification 
accuracy and many of the wetlands’ level III classes can be obtained 
if the map product is stratified into envir Graven tal units such as 
marine deposits, fluvial deposits, etc., and .analyzed by strata 
(Pectingcr, 1980). 
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Orbital remote BCUBing <iata have n level I discrimination capnb^ 
lity, ocnsionnlly allowing a level IX detailing. The information 
obtained in this work roaches some level XI details. They can be a 
first step toward many unknown Information such ns agricultural 
land use pattern and orientatiori, carrying capacity determination, 
water table level zoning, mosquito brooding cites localization 
(Cibula, 1976), and detritic sources for the estuarine food chain, 
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TABLE 1 - CLASSIFICATION SYSTEM 


1, Coal related areas 

1.1 - Deposits 

1.1.1 - Pyrite-rich coal refuse 

1.1.2 - Coal 

1.1.3 - Pyrite concentrate 
l,lv4 “ Carbonous refuse 

1.2 ~ Slurry ponds 

1.2.1 - Pyrite fine materials 

1.2.2 - Ashes 

1.2.3 - Carbonous fine materials 


2. Wetlands 

2.1 - Permanently flooded 

2.1.1 - Typha dominguensis 

2.1.2 " Spartina sp 

2.1.3 - Beach hidrosere 

2.1.4 - Stressed vegetation 

2.1.5 - floating hidrosere 

2.2 - Periodically flooded 

2.2.1 ~ Sedges plus pastures 

2.2.2 “ Pure sedges 


2.2.3 - Estiva marsh 
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3 - Agricultural lands 

3.1 - Rice 

3.2 - Pastures 

3.2.1 - Drained pastures 

3.2.2 - Frequently flooded pastures 

3.2.3 - Periodically flooded pastures 
3. 2. A “ Upland pastures 

3.3 - Bare soils 

3.3.1 “ Organic ~ humid 

3.3.2 - Organic - dry 

3.3.3 “ Sandy 
3. A Others 

A. Water bodies 
A . 1 - Lagoon 
A . 2 “ Ocean 
A. 3 - Shallow waters 

5. Urban areas 

5.1 “ Residential 

5.2 - Allotments 


5.3 - Industrial 
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6. Forestlands 

6.1 - Reforestation 

6.2 ~ Natural forests 

6.3 - Forested wetlands 

6.4 - ’’jundu'* 


^ quality 


7, Beach 

7.1 - Beaches and dunes 

7.2 - Between dunes vegetation 
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FIGURE 1 - Study area 
□ fluvial sediments 

marine sediments 

t\V| rock outcrops 



FIGURE 2 - Ground truth map 

FIGURE 3 - MSS-LANDSAT derived map 

H “ coal related areas 

0 - permanently flooded wetlands 

S - periodically flooded wetlands 
0 - vegetated agricultural lands 

+ - bare soils 

. -- water bodies 

w - urban areas 
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Figure 2 ~ Ground truth maps 
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